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1 Introduction 

In this paper we give a quantum statistical interpretation for the bracket polynomial state sum 



O 



(K) and for the Jones polynomial V/f(i). We use this quantum mechanical interpretation to give 
a new quantum algorithm for computing the Jones polynomial. This algorithm is useful for its 
conceptual simplicity and it applies to all values of the polynomial variable that lie on the unit 
circle in the complex plane. Letting C(K) denote the Hilbert space for this model, there is a 



natural unitary transformation 
such that 



U : C(K) — -> C{K) 



(K) =< i)\U\il) > 



where |?/>) is a sum over basis states for C(K). The quantum algorithm comes directly from this 
formula via the Hadamard Test. We then show that the framework for our quantum model for 
the bracket polynomial is a natural setting for Khovanov homology. The Hilbert space C(K) of 
our model has basis in one-to-one correspondence with the enhanced states of the bracket state 
summmation and is isomorphic with the chain complex for Khovanov homology with coefficients 
in the complex numbers. We show that for the Khovanov boundary operator d : C(K) — > C(K), 
we have the relationship dU + Ud = 0. Consequently, the operator U acts on the Khovanov 
homology, and we therefore obtain a direct relationship between Khovanov homology and this 



quantum algorithm for the Jones polynomial. The quantum algorithm given here is inefficient, 
and so it remains an open problem to determine better quantum algorithms that involve both the 
Jones polynomial and the Khovanov homology. 

The paper is organized as follows. Section 2 reviews the structure of the bracket polynomial, 
its state summation, the use of enhanced states and the relationship with the Jones polynomial. 
Section 3 describes the quantum statistical model for the bracket polynomial and refers to Sec- 
tion 7 (the Appendix) for the details of the Hadamard Test and the structure of the quantum 
algorithm. Section 4 describes the relationship of the quantum model with Khovanov homology 
axiomatically, without using the specific details of the Khovanov chain complex that give these 
axioms. In Section 5 we construct the Khovanov chain complex in detail and show how some of 
its properties follow uniquely from the axioms of the previous section. In Section 6 we discuss 
how the framework of this paper can be generalized to other situations where the Hilbert space 
of a quantum information system is also a chain complex. We give an example using DeRahm 
cohomology. In Section 7 (the Appendix) we detail the Hadamard test. 



2 Bracket Polynomial and Jones Polynomial 

The bracket polynomial ifTTI model for the Jones polynomial (HI HI [10l |29l is usually described 
by the expansion 

(X) = ^(X) + ^ 1 ()0 

Here the small diagrams indicate parts of otherwise identical larger knot or link diagrams. The 
two types of smoothing (local diagram with no crossing) in this formula are said to be of type A 
(A above) and type B (A^ 1 above). 

(O) = -A 2 - A- 2 
(KO) = (-A 2 -A~ 2 )(K} 
(Y) = (-A 3 )(^) 
(Y) = (-A^)(^) 

One uses these equations to normalize the invariant and make a model of the Jones polynomial. 
In the normalized version we define 

f K (A) = (-A 3 )- wr ^(K)/(0) 

where the writhe wr(K) is the sum of the oriented crossing signs for a choice of orientation of 
the link K. Since we shall not use oriented links in this paper, we refer the reader to |fTT| for 
the details about the writhe. One then has that fx (A) is invariant under the Reidemeister moves 
(again see ifTTTl ) and the original Jones polynonmial V K (t) is given by the formula 

v K (t) = Mr 1/4 ). 
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The Jones polynomial has been of great interest since its discovery in 1983 due to its relationships 
with statistical mechanics, due to its ability to often detect the difference between a knot and its 
mirror image and due to the many open problems and relationships of this invariant with other 
aspects of low dimensional topology. 

The State Summation. In order to obtain a closed formula for the bracket, we now describe it 
as a state summation. Let K be any unoriented link diagram. Define a state, S, of K to be the 
collection of planar loops resulting from a choice of smoothing for each crossing of K. There 
are two choices (A and B) for smoothing a given crossing, and thus there are 2 C ^ K > states of a 
diagram with c(K) crossings. In a state we label each smoothing with A or A^ 1 according to 
the convention indicated by the expansion formula for the bracket. These labels are the vertex 
weights of the state. There are two evaluations related to a state. The first is the product of the 
vertex weights, denoted (K\S). The second is the number of loops in the state S, denoted \ \S\\. 
Define the state summation, (K), by the formula 

(K) =J2< K\S > 5^ 

s 

where S = —A 2 — A~ 2 . This is the state expansion of the bracket. It is possible to rewrite this 
expansion in other ways. For our purposes in this paper it is more convenient to think of the loop 
evaluation as a sum of two loop evaluations, one giving —A 2 and one giving — A~ 2 . This can 
be accomplished by letting each state curve carry an extra label of +1 or —1. We describe these 
enhanced states below. 

Changing Variables. Letting c(K) denote the number of crossings in the diagram K, if we 
replace (K) by A~ c<yK ^ (K) , and then replace A 2 by — q^ 1 , the bracket is then rewritten in the 
following form: 

<x> = <x>-?< 



with (O) = (q + q" 1 )- It is useful to use this form of the bracket state sum for the sake of the 
grading in the Khovanov homology (to be described below). We shall continue to refer to the 
smoothings labeled q (or A^ 1 in the original bracket formulation) as B-smoothings. 

Using Enhanced States. We now use the convention of enhanced states where an enhanced 
state has a label of 1 or —1 on each of its component loops. We then regard the value of the loop 
q + q" 1 as the sum of the value of a circle labeled with a 1 (the value is q) added to the value of a 
circle labeled with an —1 (the value is q^ 1 ). We could have chosen the less neutral labels of +1 
and X so that 

q +1 ^+1^1 

and 

since an algebra involving 1 and X naturally appears later in relation to Khovanov homology. It 
does no harm to take this form of labeling from the beginning. 
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Consider the form of the expansion of this version of the bracket polynonmial in enhanced 
states. We have the formula as a sum over enhanced states s : 

(k) = ^(-i)^y (s) 

s 

where i(s) is the number of £>-type smoothings in s and j(s) = i(s) + A(s), with A(s) equal to 
the number of loops labeled 1 minus the number of loops labeled —1 in the enhanced state s. 

One advantage of the expression of the bracket polynomial via enhanced states is that it is 
now a sum of monomials. We shall make use of this property throughout the rest of the paper. 



3 Quantum Statistics and the Jones Polynomial 

We now use the enhanced state summation for the bracket polynomial with variable q to give a 
quantum formulation of the state sum. Let q be on the unit circle in the complex plane. (This is 
equivalent to letting the original bracket variable A be on the unit circle and equivalent to letting 
the Jones polynmial variable t be on the unit circle.) Let C(K) denote the complex vector space 
with orthonormal basis { \s) } where s runs over the enhanced states of the diagram K. The vector 
space C(K) is the (finite dimensional) Hilbert space for our quantum formulation of the Jones 
polynomial. While it is customary for a Hilbert space to be written with the letter H, we do not 
follow that convention here, due to the fact that we shall soon regard C(K) as a chain complex 
and take its homology. One can hardly avoid using % for homology. 

With q on the unit circle, we define a unitary transformation 

U : C(K) — ► C{K) 

by the formula 

U\s) = (-l)^V' (s) |s) 

for each enhanced state s. Here i(s) and j(s) are as defined in the previous section of this paper. 
Let 

s 

The state vector \ip) is the sum over the basis states of our Hilbert space C(K). For convenience, 
we do not normalize to length one in the Hilbert space C(K). We then have the 

Lemma. The evaluation of the bracket polynomial is given by the following formula 

(K) = M W>. 

Proof. 

MW> =EE( s 'l(- 1 ) ,:(s V (s) k) = EE(- 1 ) !(S) 9 J(S) (^) 

s' s s i s 
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= E(- i r ( v (s) = (io, 

s 

since 

(s'\s) = 6(s,s') 

where S(s, s') is the Kronecker delta, equal to 1 when s = s' and equal to otherwise. // 

Thus the bracket polyomial evaluation is a quantum amplitude for the measurement of the 
state U\ip) in the (tp\ direction. Since (ip\U\ip) can be regarded as a diagonal element of the 
transformation U with respect to a basis containing this formula can be taken as the foun- 
dation for a quantum algorithm that computes the bracket of K via the Hadamard test. See the 
Appendix to this paper, for a discussion of the Hadamard test and the corresponding quantum 
algorithm. 

A few words about quantum algorithms are appropropiate at this point. In general one begins 
with an intial state \ifj) and a unitary transformation U. Quantum processes are modeled by uni- 
tary transformations, and it is in principle possible to create a physical process corresponding to 
any given unitary transformation. In practice, one is limited by the dimensions of the spaces in- 
volved and by the fact that physical quantum states are very delicate and subject to decoherence. 
Nevertheless, one designs quantum algorithms at first by finding unitary operators that represent 
the informaton in the problem one wishes to calculate. Then the quantum part of the quantum 
computation is the physical process that produces the state U\ip) from the initial state . At this 
point U\ijj) is available for measurement. Measurement means interaction with the enviroment, 
and will happen in any case, but one intends a controlled circumstance in which the measurement 
can occur. If {|ei), • ■ ■ \e n )} is a basis for the Hilbert space (here denoted by C(K)) then we 
would have 

U\i>) = H^fc|efc), 

k 

a linear combination of the basis elements. On measuring with respect to this basis, each \e^) cor- 
responds to an observable outcome and this outcome will occur with frequency \zk\ 2 /(Y,i \ z i\ 2 )- 
The form of quantum computation consists in repeatedly preparing and running the process to 
find the frequency with which certain key observations occur. In the case of our algorithm we 
are interested in the frequency of measuring \tp) itself and this corresponds to the inner product 
(ip\U\ip). However, a direct measurement from U\tp) will yield only approximations to the abso- 
lute square of (ip\U\iJ)). For this reason we use the more refined Hadamard test as described in the 
Appendix. There is much more to be said about quantum algorithms in general and about quan- 
tum algorithms for the Jones polynomial in particular. The reader can examine lfTl [T4l[T5l[T6ll26l 
for more information. 

It is useful to formalize the bracket evaluation as a quantum amplitude. This is a direct way 
to give a physical interpretation of the bracket state sum and the Jones polynomial. Just how 
this process can be implemented physically depends upon the interpretation of the Hilbert space 
C(K). It is common practice in theorizing about quantum computing and quantum information 
to define a Hilbert space in terms of some mathematically convenient basis (such as the enhanced 
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states of the knot or link diagram K ) and leave open the possibility of a realization of the space 
and the quantum evolution operators that have been defined upon it. In principle any finite dimen- 
sional unitary operator can be realized by some physical system. In practice, this is the problem of 
constructing quantum computers and quantum information devices. It is not so easy to construct 
what can be done in principle, and the quantum states that are produced may be all too short-lived 
to produce reliable computation. Nevertheless, one has the freedom to create spaces and opera- 
tors on the mathematical level and to conceptualize these in a quantum mechanical framework. 
The resulting structures may be realized in nature and in present or future technology. In the case 
of our Hilbert space associated with the bracket state sum and its corresponding unitary trans- 
formation U, there is rich extra structure related to Khovanov homology that we discuss in the 
next section. One hopes that in a (future) realization of these spaces and operators, the Khovanov 
homology will play a key role in quantum information related to the knot or link K. 

There are a number of conclusions that we can draw from the formula (K) = (tp\U\tp). 
First of all, this formulation constitutes a quantum algorithm for the computation of the bracket 
polynomial (and hence the Jones polynomial) at any specialization where the variable is on the 
unit circle. We have defined a unitary transformation U and then shown that the bracket is an 
evaluation in the form (ip\U\iJj). This evaluation can be computed via the Hadamard test E4ll 
and this gives the desired quantum algorithm. Once the unitary transformation is given as a 
physical construction, the algorithm will be as efficient as any application of the Hadamard test. 
The present algorithm requires an exponentially increasing complexity of construction for the 
associated unitary transformation, since the dimension of the Hilbert space is equal to the 2 e ^ 
where e(K) is the number of enhanced states of the diagram K. (Note that e(K) = Y,s^ S ^ 
where S runs over the 2 C ( K > standard bracket states, c(K) is the number of crossings in the 
diagram and \ \S\ \ is the number of loops in the state S.) Nevertheless, it is significant that the 
Jones polynomial can be formulated in such a direct way in terms of a quantum algorithm. By 
the same token, we can take the basic result of Khovanov homology that says that the bracket is 
a graded Euler characteristic of the Khovanov homology as telling us that we are taking a step in 
the direction of a quantum algorithm for the Khovanov homology itself. This will be discussed 
below. 



4 Khovanov Homology and a Quantum Model for the Jones 
Polynomial 

In this section we outline how the Khovanov homology is related with our quantum model. This 
can be done essentially axiomatically, without giving the details of the Khovanov construction. 
We give these details in the next section. The outline is as follows: 

1 . There is a boundary operator d defined on the Hilbert space of enhanced states of a link 
diagram K 

d : C(K) — > C(K) 
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such that dd = and so that if C M = C l ' j (K) denotes the subspace of C(K) spanned by 
enhanced states \s) with i — i(s) and j = j(s), then 



d : C j — > C i+1 ' j . 

That is, we have the formulas 

#)) = i(l«)) + l 

and 

for each enhanced state s. In the next section, we shall explain how the boundary operator 
is constructed. 

2. Lemma. By defining U : C(K) — > C(K) as in the previous section, via 

U\s) = (-l) i(s V (s) |s>, 
we have the following basic relationship between U and the boundary operator d : 

Ud + dU = 0. 

Proof. This follows at once from the definition of U and the fact that d preserves j and 
increases % to % + 1. // 

3. From this Lemma we conclude that the operator U acts on the homology of C(K). We 
can regard H(C(K)) = Ker(d) / Image(d) as a new Hilbert space on which the unitary 
operator U acts. In this way, the Khovanov homology and its relationship with the Jones 
polynomial has a natural quantum context. 

4. For a fixed value of j, 

C'- J = ®iC ij 

is a subcomplex of C(K) with the boundary operator d. Consequently, we can speak of 
the homology iJ(C* J ). Note that the dimension of C* J is equal to the number of enhanced 
states \s) with i = i(s) and j = j(s). Consequently, we have 

(k) = £^)(-i)*o = 

s j i 

j j 

Here we use the definition of the Euler characteristic of a chain complex 

x (C'j) = J2(-iydirn(C ij ) 



and the fact that the Euler characteristic of the complex is equal to the Euler characteristic 
of its homology. The quantum amplitude associated with the operator U is given in terms 
of the Euler characteristics of the Khovanov homology of the link K. 

(K) = ty\U\1>) =J2v j x(H(C'J(K))). 

j 

Our reformulation of the bracket polynomial in terms of the unitary operator U leads to a 
new viewpoint on the Khovanov homology as a representation space for the action of U. The 
bracket polynomial is then a quantum amplitude that expressing the Euler characteristics of the 
homology associated with this action. The decomposition of the chain complex into the parts 
C l 'i(K) corresponds to the eigenspace decomposition of the operator U. The reader will note 
that in this case the operator U is already diagonal in the basis of enhanced states for the chain 
complex C(K). We regard this reformulation as a guide to further questions about the relationship 
of the Khovanov homology with quantum information associated with the link K. 

As we shall see in the next section, the internal combinatorial structure of the set of enhanced 
states for the bracket summation leads to the Khovanov homology theory, whose graded Euler 
characteristic yields the bracket state sum. Thus we have a quantum statistical interpretation of 
the Euler characteristics of the Khovanov homology theory, and a conceptual puzzle about the 
nature of this relationship with the Hilbert space of that quantum theory. It is that relationship 
that is the subject of this paper. The unusual point about the Hilbert space is the each of its basis 
elements has a specific combinatorial structure that is related to the topology of the knot K . Thus 
this Hilbert space is, from the point of view of its basis elements, a form of "taking apart" of the 
topological structure of the knot that we are interested in studying. 

Homological structure of the unitary transformation. We now prove a general result about 
the structure of a chain complex that is also a finite dimensional Hilbert space. Let C be a chain 
complex over the complex numbers with boundary operator 

d : C — ■> C i+ \ 

with C denoting the direct sum of all the C\ % — 0, 1, 2, • • • n (for some n). Let 

U:C^C 

be a unitary operator that satisfies the equation U d+dU = 0. We do not assume a second grading 
j as occurs in the Khovanov homology. However, since U is unitary, it follows [|22ll that there is a 
basis for C in which U is diagonal. Let B = {\s)} denote this basis. Let A s denote the eigenvalue 
of U corresponding to \s) so that U\s) = X s \s). Let a S)S / be the matrix element for d so that 

d\s) = J2 a s,s'\s') 

s' 

where s' runs over a set of basis elements so that i(s') = i(s) + 1. 
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Lemma. With the above conventions, we have that for \s') a basis element such that a S)S > ^ 
then A s / = — \ s . 
Proof. Note that 

Ud\s) = U(J2 a.y |«'))=E A *' I s '> 

s' s' 

while 

dU\s) = d\ s \s) = ^a SjS /A s |s'). 

s' 

Since U d+dU = 0, the conclusion of the Lemma follows from the independence of the elements 
in the basis for the Hilbert space. // 

In this way we see that eigenvalues will propagate forward from C° with alternating signs ac- 
cording to the appearance of successive basis elements in the boundary formulas for the chain 
complex. Various states of affairs are possible in general, with new eigenvaluues starting at some 
C k for k > 0. The simplest state of affairs would be if all the possible eigenvalues (up to multi- 
plication by —1) for U occurred in C° so that 

where A runs over all the distinct eigenvalues of U restricted to C°, and C° is spanned by all 
\s) in C° with U\s) = X\s). Let us take the further assumption that for each A as above, the 
subcomplexes 

A • U A f U -A f U +A f ' ' ' L -(-l)"A 

have C = ©aC* as their direct sum. With this assumption about the chain complex, define 
\ip) = J2 S I s ) as before, with \s) running over the whole basis for C. Then it follows just as in the 
beginning of this section that 

<VW> = E*xW?x))- 

A 

Here \ denotes the Euler characteristic of the homology. The point is, that this formula for 
(ip\U\ip) takes exactly the form we had for the special case of Khovanov homology (with eigen- 
values (— 1) V ), but here the formula occurs just in terms of the eigenspace decomposition of the 
unitary transformation U in relation to the chain complex. Clearly there is more work to be done 
here and we will return to it in a subsequent paper. 

Remark on the density matrix. Given the state [ip), we can define the density matrix 

P=hA>(# 

With this definition it is immediate that 

Tr{Up) = (i/j\U\ip) 

where Tr(M) denotes the trace of a matrix M. Thus we can restate the form of our result about 
Euler characteristics as 

Tr(Up)=J2 X x(H(Cl))- 

A 
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In searching for an interpretation of the Khovanov complex in this quantum context it is useful to 
use this reformulation. For the bracket we have 

(K) =< ifj\U\ij >= Tr{Up). 

Remark on physics. We are taking the whole Hilbert space as the state space of a physical sys- 
tem. If the system were a classical system then the energetic states of the classical system would 
correspond to the basis we have chosen for the Hilbert space. This is in line with the concepts of 
quantum mechanics, since the classical states are then in correspondence with a basis for mea- 
surement. In quantum context we think of the elements in the Hilbert space as corresponding 
to superpositions of possible measurements. All this is then transposed to topological configu- 
rations for the knot or link. But the new ingredient is the one from topology that will make the 
states of the knot into a chain complex and lead to homological computations. Since we are using 
enhanced states, each state is a generator for the chain complex. Thus we can regard C(K) itself 
as the chain complex (over the complex numbers) and add in extra grading structure to boot. 
There is also the fine structure of the underlying state circles, but this is not seen by the Hilbert 
space itself and the chain complex structure can certainly be written just in terms of the basis of 
the Hilbert space. We want to know what is the relationship between the unitary transformation 
and the homological structure. This relation is given by the way the grading works in relation to 
(— lyqi for each state: the constancy of the quantum grading j under the action of the differen- 
tial. It is this that makes (il>\U\i(j) a graded Euler characteristic for the homology. In this way the 
unitary transformation is linked with the structure of state transitions that govern the homology. 
At the end of this section we indicate a more general approach to this pattern. 

The key conceptual issue is the construction of a Hilbert space whose basis is the set of 
observable states of a physical system. We can always do this for a statistical mechanical system. 
Let s denote such states and E(s) denote the energy of a state s (observable and real). Define a 
unitary transformation U (t) (where t denotes time) by 

U(t)\s) = e {it/n)E{s) \s). 

Define the amplitude 

A(t) = E s e {lt/n)E{s) . 

Note that we have the formula 

A(t) = (^\u{tm ) 

where \ip) = E s |s) just as before. Now U(t) can be regarded as the quantum evolution of the 
state of the system and we have written this amplitude in analogy with a partition function for the 
statistical mechanics model. 

One more conceptual problem: We consider an amplitude of the form 

y,m) = x;(-i)*v (s) 
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where q = e is a unit complex number. The question is, how does this fit the pattern 

A(t) = Y / e {u/h)E{s) ? 

s 

So we want 

(t/h)E(s) = m(s) + 9j(s) 

which gives 

E(s) = {h/t)(m(8) + 0j( 8 )). 

Evolution from time t = is out of the question. Probably the best interpretation is to think of 
the evolution starting at t — 1 . 

In this point of view, the Hilbert space for expressing the bracket polynomial as a quantum 
statistical amplitude is quite naturally the chain complex for Khovanov homology with complex 
coefficients, and the unitary transformation that is the structure of the bracket polynomial acts on 
the homology of this chain complex. This means that the homology classes contain information 
preserved by the quantum process that underlies the bracket polynomial. We would like to exploit 
this direct relationship between the quantum model and the Khovanov homology to obtain deeper 
information about the relationship of topology and quantum information theory, and we would 
like to use this relationship to probe the properties of these topological invariants. 



5 Background on Khovanov Homology 

In this section, we describe Khovanov homology along the lines of Ifl9l l2l. and we tell the story 
so that the gradings and the structure of the differential emerge in a natural way. This approach to 
motivating the Khovanov homology uses elements of Khovanov's original approach, Viro's use of 
enhanced states for the bracket polynomial [|28l . and Bar-Natan's emphasis on tangle cobordisms 
||3l . We use similar considerations in our paper IfTTll . 

Two key motivating ideas are involved in finding the Khovanov invariant. First of all, one 
would like to categorify a link polynomial such as (K) . There are many meanings to the term 
categorify, but here the quest is to find a way to express the link polynomial as a graded Euler 
characteristic (K) = Xq{H(K)) for some homology theory associated with (K). 

We will use the bracket polynomial and its enhanced states as described in the previous sec- 
tions of this paper. To see how the Khovanov grading arises, consider the form of the expansion 
of this version of the bracket polynomial in enhanced states. We have the formula as a sum over 
enhanced states s : 

(K) = ^(-l) i(s V (s) 
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where i(s) is the number of £>-type smoothings in s, A(s) is the number of loops in s labeled 
1 minus the number of loops labeled X, and j(s) = i(s) + A(s). This can be rewritten in the 
following form: 

(K) = Y,(-iy<l J dim(C ij ) 

i,j 

where we define C JJ to be the linear span (over the complex numbers for the purpose of this paper, 
but over the integers or the integers modulo two for other contexts) of the set of enhanced states 
with i(s) = i and j(s) = j. Then the number of such states is the dimension dim{C^). 

We would like to have a bigraded complex composed of the C y with a differential 

d : C lj — ► C i+lj . 

The differential should increase the homological grading i by 1 and preserve the quantum grading 
j. Then we could write 

(k) = £^'E(-i)^™(n =xyx(c" i ), 

3 i 3 

where x{C^) is the Euler characteristic of the subcomplex C ,J for a fixed value of j. 

This formula would constitute a categorification of the bracket polynomial. Below, we shall see 
how the original Khovanov differential d is uniquely determined by the restriction that j(ds) = 
j(s) for each enhanced state s. Since j is preserved by the differential, these subcomplexes C' j 
have their own Euler characteristics and homology. We have 

where HiC* '-*) denotes the homology of the complex C*K We can write 

j 

The last formula expresses the bracket polynomial as a graded Euler characteristic of a homology 
theory associated with the enhanced states of the bracket state summation. This is the categorifi- 
cation of the bracket polynomial. Khovanov proves that this homology theory is an invariant of 
knots and links (via the Reidemeister moves of Figure 1), creating a new and stronger invariant 
than the original Jones polynomial. 

We will construct the differential in this complex first for mod-2 coefficients. The differential 
is based on regarding two states as adjacent if one differs from the other by a single smoothing 
at some site. Thus if (s, r) denotes a pair consisting in an enhanced state s and site r of that 
state with r of type A, then we consider all enhanced states s' obtained from s by smoothing at 
r and relabeling only those loops that are affected by the resmoothing. Call this set of enhanced 
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Figure 1 : Reidemeister Moves 



states S'[s, r]. Then we shall define the partial differential d T (s) as a sum over certain elements 
in S'[s, r], and the differential by the formula 



with the sum over all type A sites r in s. It then remains to see what are the possibilities for d T (s) 
so that j(s) is preserved. 

Note that if s' E S'[s, r], then i(s') = i(s) + 1. Thus 

From this we conclude that j(s) = j(s') if and only if A(s') = A(s) — 1. Recall that 



where [s : +] is the number of loops in s labeled +1, [s : — ] is the number of loops labeled — 1 
(same as labeled with X) and j(s) = i(s) + A(s). 

Proposition. The partial differentials d T (s) are uniquely determined by the condition that j(s') = 
j(s) for all s' involved in the action of the partial differential on the enhanced state s. This unique 
form of the partial differential can be described by the following structures of multiplication and 
comultiplication on the algebra A = k[X]/(X 2 ) where k = Z/2Z for mod-2 coefficients, or 
k = Z for integral coefficients. 

1. The element 1 is a multiplicative unit and X 2 = 0. 

2. A(l) = 1 ® X + X g> 1 and A(X) = X®X. 

These rules describe the local relabeling process for loops in a state. Multiplication corresponds 
to the case where two loops merge to a single loop, while comultiplication corresponds to the 
case where one loop bifurcates into two loops. 



d{s)=Y,dr{s) 



T 



\{s) = [s:+]-[s:-} 
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Proof. Using the above description of the differential, suppose that there are two loops at r that 
merge in the smoothing. If both loops are labeled 1 in s then the local contribution to A(s) is 2. 
Let s' denote a smoothing in S[s, r]. In order for the local A contribution to become 1, we see 
that the merged loop must be labeled 1. Similarly if the two loops are labeled 1 and X, then the 
merged loop must be labeled X so that the local contribution for A goes from to —1. Finally, 
if the two loops are labeled X and X, then there is no label available for a single loop that will 
give —3, so we define d to be zero in this case. We can summarize the result by saying that there 
is a multiplicative structure m such that m(l, 1) = l,m(l,X) = m(X, 1) = x, m(X, X) = 0, 
and this multiplication describes the structure of the partial differential when two loops merge. 
Since this is the multiplicative structure of the algebra A = k[X]/ (X 2 ), we take this algebra as 
summarizing the differential. 

Now consider the case where s has a single loop at the site r. Smoothing produces two loops. 
If the single loop is labeled X, then we must label each of the two loops by X in order to make A 
decrease by 1. If the single loop is labeled 1, then we can label the two loops by X and 1 in either 
order. In this second case we take the partial differential of s to be the sum of these two labeled 
states. This structure can be described by taking a coproduct structure with A(X) = X <g> X and 
A(l) = l(g)X+X<g)l. We now have the algebra A = k[X]/(X 2 ) with product m : A® A — > A 
and coproduct A : A — ► A® A, describing the differential completely. This completes the 
proof. // 

Partial differentials are defined on each enhanced state s and a site r of type A in that state. 
We consider states obtained from the given state by smoothing the given site r. The result of 
smoothing r is to produce a new state s' with one more site of type B than s. Forming s' from s 
we either amalgamate two loops to a single loop at r, or we divide a loop at r into two distinct 
loops. In the case of amalgamation, the new state s acquires the label on the amalgamated circle 
that is the product of the labels on the two circles that are its ancestors in s. This case of the 
partial differential is described by the multiplication in the algebra. If one circle becomes two 
circles, then we apply the coproduct. Thus if the circle is labeled X, then the resultant two 
circles are each labeled X corresponding to A(X) = X <g> X. If the orginal circle is labeled 
1 then we take the partial boundary to be a sum of two enhanced states with labels 1 and X 
in one case, and labels X and 1 in the other case, on the respective circles. This corresponds to 
A(l) = 1®X+X®1. Modulo two, the boundary of an enhanced state is the sum, overall sites of 
type A in the state, of the partial boundaries at these sites. It is not hard to verify directly that the 
square of the boundary mapping is zero (this is the identity of mixed partials!) and that it behaves 
as advertised, keeping j(s) constant. There is more to say about the nature of this construction 
with respect to Frobenius algebras and tangle cobordisms. In Figures 2,3 and 4 we illustrate how 
the partial boundaries can be conceptualized in terms of surface cobordisms. The equality of 
mixed partials corresponds to topological equivalence of the corresponding surface cobordisms, 
and to the relationships between Frobenius algebras and the surface cobordism category. In 
particular, in Figure 4 we show how in a key case of two sites (labeled 1 and 2 in that Figure) the 
two orders of partial boundary are 

d 2 di = (1 <g> m) o (A ® 1) 
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Figure 2: SaddlePoints and State Smoothings 

and 

8182 = A o m. 

In the Frobenius algebra A = k[X]/ (X 2 ) we have the identity 

(1 ® m) o (A (8) 1) = A om. 

Thus the Frobenius algebra implies the identity of the mixed partials. Furthermore, in Figure 
3 we see that this identity corresponds to the topological equivalence of cobordisms under an 
exchange of saddle points. There is more to say about all of this, but we will stop here. The 
proof of invariance of Khovanov homology with respect to the Reidemeister moves (respecting 
grading changes) will not be given here. See lfT9l l2ll3l. It is remarkable that this version of Kho- 
vanov homology is uniquely specified by natural ideas about adjacency of states in the bracket 
polynomial. 

Remark on Integral Differentials. Choose an ordering for the crossings in the link diagram K 
and denote them by 1, 2, ■ • • n. Let s be any enhanced state of K and let 8i(s) denote the chain 
obtained from s by applying a partial boundary at the i-th site of s. If the i-th site is a smoothing 
of type A^ 1 , then 8i(s) — 0. If the i-th site is a smoothing of type A, then <9;(s) is given by the 
rules discussed above (with the same signs). The compatibility conditions that we have discussed 
show that partials commute in the sense that 8i (8j(s) ) = 8j(di(s)) for all i and j. One then defines 
signed boundary formulas in the usual way of algebraic topology. One way to think of this regards 
the complex as the analogue of a complex in de Rahm cohomology. Let {dxi, dx 2 , ■ ■ • , dx n } be 
a formal basis for a Grassmann algebra so that dxi A dxj = —dxj A dx,i Starting with enhanced 
states s in C°(K) (that is, state with all A-type smoothings) Define formally, di(s) = 8i{s)dxi 
and regard di(s) as identical with 8i(s) as we have previously regarded it in C l {K). In general, 
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F G 

Figure 3: Surface Cobordisms 




9 2 9i = d(g)m)(A(g)1) 
9i9 2 = (A)(m) 

9i9 2 =9 2 9i 



Figure 4: Local Boundaries Commute 
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given an enhanced state s in C k (K) with 5-smoothings at locations z'i < i 2 < ■ ■ ■ < ik, we 
represent this chain as s dx ix A ■ ■ ■ A dx ik and define 

n 

d(s dxi x A ■ ■ ■ A dx ik ) = <9?(s) dxj A da^ A ■ ■ ■ A o?Xi fe , 

i=i 

just as in a de Rahm complex. The Grassmann algebra automatically computes the correct signs 
in the chain complex, and this boundary formula gives the original boundary formula when we 
take coefficients modulo two. Note, that in this formalism, partial differentials <9« of enhanced 
states with a £>-smoothing at the site i are zero due to the fact that dx. t A dxi = in the Grass- 
mann algebra. There is more to discuss about the use of Grassmann algebra in this context. For 
example, this approach clarifies parts of the construction in lfT8l . 

It of interest to examine this analogy between the Khovanov (co)homology and de Rahm 
cohomology. In that analogy the enhanced states correspond to the differentiable functions on a 
manifold. The Khovanov complex C k (K) is generated by elements of the form s dx. lx A • • -Adxi k 
where the enhanced state s has 5-smoothings at exactly the sites z'i, • • • , z*.. If we were to follow 
the analogy with de Rahm cohomology literally, we would define a new complex DR(K) where 
DR k (K) is generated by elements s dx^ A • • • A dx ik where s is any enhanced state of the link 
K. The partial boundaries are defined in the same way as before and the global boundary formula 
is just as we have written it above. This gives a new chain complex associated with the link K. 
Whether its homology contains new topological information about the link K will be the subject 
of a subsequent paper. 

A further remark on de Rham cohomology. There is another deep relation with the de Rham 
complex: In j]2Tfl it was observed that Khovanov homology is related to Hochschild homology 
and Hochschild homology is thought to be an algebraic version of de Rham chain complex (cyclic 
cohomology corresponds to de Rham cohomology), compare ll23l . 



6 Other Homological States 

The formalism that we have pursued in this paper to relate Khovanov homology and the Jones 
polynomial with quantum statistics can be generalized to apply to other situations. It is possible to 
associate a finite dimensional Hilbert space that is also a chain complex (or cochain complex) to 
topological structures other than knots and links. To formulate this, let X denote the topological 
structure and C(X) denote the associated linear space endowed with a boundary operator d : 
C(X) — > C(X), with dd = 0. We want to consider situations where there is a unitary operator 
U : C(X) — > C(X) such that Ud + dU = so that U induces a unitary action on T-L(X), the 
homology of C(X) with respect to d. 



17 



For example, let M be a differentiable manifold and C(M) denote the DeRham complex of 
M over the complex numbers. Then for a differential form of the type f(x)u in local coordinates 
Xi, ■ ■ ■ , x n and cu a wedge product of a subset of dx\ ■ ■ ■ dx n , we have 

n 

d(fu) = J2( d f /dxi)dxi A to. 

i=l 

Here d is the differential for the DeRham complex. Then C(M) has basis the set of \f(x)co) 
where to = dx, lx A ■ • • A dxi k with i\ < ■ ■ ■ < iy.. We could achieve Ud + dll = if U is a 
very simple unitary operator (e.g. multiplication by phases that do not depend on the coordinates 
Xi) but in general it will be an interesting problem to determine all unitary operators U with this 
property. 

Even in the case of Khovanov homology, we could keep the homology theory fixed and ask 
for other unitary operators U that satisfy Ud + dU = 0. Knowing other examples of such oper- 
ators would shed light on the nature of Khovanov homology from the point of view of quantum 
statistics. 



7 Appendix - The Hadamard Test 

In order to make a quantum computation of the trace of a unitary matrix U, one can use the 
Hadamard test to obtain the diagonal matrix elements (ip\U\ijj) of U. The trace is then the sum of 
these matrix elements as runs over an orthonormal basis for the vector space. In the applica- 
tion to the algorithm described here for the Jones polynomial it is only necessary to compute one 
number of the form (^\U\^). The Hadamard test proceeds as follows. 

We first obtain 

as an expectation by applying the Hadamard gate H 

#|0} = ^(|0> + |1}) 



to the first qubit of 



#|1) = ^(|0)-|1» 



Cu o (H ® 1)|0M = -J=(|0) ® M + |1> ® W>- 



Here Cu denotes controlled U, acting as U when the control bit is |1) and the identity mapping 
when the control bit is |0). We measure the expectation for the first qubit |0) of the resulting state 

(H ® 1) o Cu o (H ® l)\0)\i/>) = 
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Figure 5: Hadamard Test 



|c^io> ® m +h\i)® um = |((io> + |i» ® i^) + do) - |i» ® 



This expectation is 



2(|0)®(|V) + W» + |l)®(lV'>- 



In Figure 5, we illustrate this computation with a diagram that indicates the structure of the 
test with parallel lines corresponding to tensor products of the single qubit space (with three lines 
chosen for illustration as the size of U). The extra tensor factor is indicated on the top line with the 
Hadamard matrix H indicated by a box and the control of U indicated by a circle with a vertical 
control line extending down to the U -box. The half -circle on the top line on the right stands for 
the measurement of that line that is used for the computation. Thus Figure 5 represents a circuit 
diagram for the quantum computation of \ + \Re(ip\U\il)) and hence the quantum computation 
of | + ^Re(K) when U is taken to be the unitary tranformation corresponding to the bracket 
polynomial, as discussed in previous sections of this paper. 



The imaginary part is obtained by applying the same procedure to 

^=(|0>® |V)-i|l>® W> 

This is the method used in 03 Q31 S3, and the reader may wish to contemplate its efficiency in 
the context of this simple model. Note that the Hadamard test enables this quantum computation 
to estimate the trace of any unitary matrix U by repeated trials that estimate individual matrix 
entries (ip\U\ip). 



References 

[1] D. Aharonov, V. Jones, Z. Landau, A polynomial quantum algorithm for approximating the 
Jones polynomial, quant-ph/0511096, 



19 



[2] Bar-Natan, D. (2002), On Khovanov's categorification of the Jones polynomial, Algebraic 
and Geometric Topology, 2(16), pp. 337-370. 

[3] Bar-Natan, D. (2005), Khovanov's homology for tangles and cobordisms, Geometry and 
Topology, 9-33, pp. 1465-1499. arXiv:mat .GT/0410495 

[4] R.J. Baxter. Exactly Solved Models in Statistical Mechanics. Acad. Press (1982). 

[5] R. Dijkgraaf, D. Orland and S. Reffert, Dimer models, free fermions and super quantum 
mechanics, larXiv : 0705 . 1 6451 

[6] S. Gukov, Surface operators and knot homologies, larXiv: 0706 . 23691 

[7] L.Helme-Guizon, J.H.Przytycki, Y.Rong, Torsion in Graph Homology, Fundamenta Math- 
ematicae, 190; June 2006, 139-177. 

[8] V.F.R. Jones, A polynomial invariant for links via von Neumann algebras, Bull. Amer. Math. 
Soc. 129(1985), 103-112. 

[9] V.F.R.Jones. Hecke algebra representations of braid groups and link polynomials. Ann. of 
Math. 126 (1987), pp. 335-338. 

[10] V.F.R.Jones. On knot invariants related to some statistical mechanics models. Pacific J. 
Math., vol. 137, no. 2 (1989), pp. 311-334. 

[11] L.H. Kauffman, State models and the Jones polynomial, Topology 26 (1987), 395-407. 

[12] L.H. Kauffman, Statistical mechanics and the Jones polynomial, AMS Contemp. Math. 
Series 78 (1989), 263-297. 

[13] L.H. Kauffman, Knots and Physics, World Scientific Publishers (1991), Second Edition 
(1993), Third Edition (2002). 

[14] L.H. Kauffman, Quantum computing and the Jones polynomial, in Quantum Computa- 
tion and Information, S. Lomonaco, Jr. (ed.), AMS CONM/305, 2002, pp. 101-137. 

math. QA/ 01052551 

[15] L. H. Kauffman and S. Lomonaco Jr., A Three-stranded quantum algorithm for the Jones 
polynonmial, in "Quantum Information and Quantum Computation V", Proceedings of 
Spie, April 2007, edited by E.J. Donkor, A.R. Pirich and H.E. Brandt, pp. 65730T1-17, 
Intl Soc. Opt. Eng. 

[16] S. J. Lomonaco Jr. and L. H. Kauffman, Quantum Knots and Mosaics, Jour- 
nal of Quantum Information Processing, Vol. 7, Nos. 2-3, (2008), pp. 85 - 115. 

arxiv.org/abs/0 8 05. 033 9 

[17] H.A. Dye, L.H. Kauffman, V.O. Manturov, On two categorifications of the arrow polyno- 
mial for virtual knots, larXiv: 0906 . 34081 



20 



[18] V.O. Manturov, Khovanov homology for virtual links with arbitrary coefficients, 

|math.GT/0 601l"52l 

[19] Khovanov, M. (1997), A categorification of the Jones polynomial, Duke Math. 7,101 (3), 
pp.359-426. 

[20] Lee, T. D., Yang, C. N. (1952), Statistical theory of equations of state and phase transitions 
II., lattice gas and Ising model, Physical Review Letters, 87, 410-419. 

[21] J. H. Przytycki, When the theories meet: Khovanov homology as Hochschild homology of 
links, larXiv: math. GT/ 50 9334| 

[22] S. Lang, Algebra, Springer- Verlag, New York, Inc. (2002). 

[23] J-L. Loday, Cyclic Homology, Grund. Math. Wissen. Band 301, Springer- Verlag, Berlin, 
1992 (second edition, 1998). 

[24] M. A. Nielsen and I. L. Chuang, "Quantum Computation and Quantum Information", Cam- 
bridge University Press (2000). 

[25] E.F. Jasso-Hernandez and Y. Rong, A categorification of the Tutte polynomial, Alg. Geom. 
Topology 6 2006, 2031-2049. 

[26] Raimund Marx, Amr Fahmy, Louis Kauffman, Samuel Lomonaco, Andreas Sporl, Nikolas 
Pomplun, John Myers, Steffen J. Glaser, NMR Quantum Calculations of the Jones Polyno- 
mial, larXiv: 0909 . 10801 

[27] M. Stosic, Categorification of the dichromatic polynomial for graphs, 
larXiv:math/0504239l v2. JKTR, Vol. 17, No. 1, (2008), pp. 31-45. 

[28] O. Viro (2004), Khovanov homology, its definitions and ramifications, Fund. Math., 184 
(2004), pp. 317-342. 

[29] E. Witten. Quantum Field Theory and the Jones Polynomial. Comm. in Math. Phys. Vol. 
121 (1989), 351-399. 



21 



